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Abstract 

Traditional large eddy simulation is based on Kolmogrov's hypothesis, and done in the inertial range. In 
incrtial range the LES model coefficient is scale- invariant. In many cases, such as computing in the boundary 
layer, the filter scale is not in the inertial range, and the model coefficient changes drastically with the filter 
scale, thus the assumption of scale-invariance is not proper. In this paper, we propose a general principle of 
constraint, and then we deduce two concrete constraints based on the existing spectrum theory and priori 
data. The constrained dynamic Smogorinsky model can predict energy spectrum well in different range of 
forced and decayed turbulence. Dynamic mixed models with one of the constraints shows more efficient 
than some classical models in predicting energy spectrum and SGS dissipation, and also displays a strong 
correlation with the real stress in different range of the turbulent flows. 
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1 



I. INTRODUCTION 



Large eddy simulation (LES) is an important method in numerical simulation of turbulent 
flows pHa]. Traditional large eddy simulation is based on the K41 theory (Kolmogrov,1941), and 
done in the inertial range. In the inertial range, the turbulent motions are statistically isotropic, 
and the statistics of the motions are decided by the dissipation rate and wave number .When 
doing LES, the motions larger than filter scale are solved directly, and the motions smaller than 
filter scale (SGS) are modeled by the resolved quantities. Several kinds of SGS models have been 
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Recently, a 



proposed so far, such as eddy-viscosity model, dynamic model, vortex model, etc. 
type of constrained LES method has already been developed Q|. The traditional S 
the assumption of scale- invariance in the inertial range. Dynamic Smogorinsky model (DSM)is a 
classic model in large eddy simulation of turbulent flows 
the Germano identity, 



. Here dynamic procedure is based on 



L{j — T{j T{j — V>i Uj Ui Uj , ( 1 ) 

Where Ly is called the resolved stress, Tij = UjUj — UiUj is the SGS stress at the filtered scale 



A, and = UjUj — UiUj is the SGS stress at the filtered scale aA(a > 1). Here a tilde denotes 
filtering at scale A and a bar denotes test-filtering at scale at a A. In ([T]), taking Smogorinsky 
model stress r™ od and T™ d to replace Tij and T^, we can get 



^mod _ Juried _ T mod 



Where LYT A is called the resolved stress, 



T™« = C sm A'\S\S tJ , (3) 

J^ od =C sm (aA) 2 \^S ij , (4) 

~ . / ~ ~ \ X /2 

and here Sij is resolved strain tensor at filtering scale A, \S\ = {2SijSij\ . A square error 
between Ly and L™ od can be given by 

F n°d =(( ^._ L mod ) 2 )) (5) 

Here the bracket (•) denotes ensemble average of the turbulence flow. According to the hypothesis 



of scale-invariance in the inertial range of turbulence flow, the model coefficients are same in t vj 



mod 

3 ai c odinc 111 / , 

and TP? od . Minimizing ([5]), we can get the expression of the model coefficients, 



m ~ (MijMij) ' W 



where, 



Mij = (aA) 2 ^ - A 2 |5|^, 



(7) 



When do LES using the dynamic Smogorinsky model, the test-scale and even the filter-scale are 



often out of the inertial range. Special 



y, more and more researches pay attention to the behavior 



of turbulence in dissipation range jij-llll]. Thus the assumption of scale- invariance is not proper yet 
[12I ]. and the validity of the traditional LES result is also not credible. In order to solve this problem, 
several methods have been proposed so far. Ronchi et al. have proposed a complex dynamic model 
[13] ]. and the model coefficient varies with scale. Although the practicability of the model has not 
been tested widely, it supplies a feasible idea to us all. Using several existing energy spectrum, 
P. R. Voke gives the approximative coefficient of Smogorinsky model in dissipation range [rJ. 
Recently, Meneveau et al. have proposed a scale-dependent dynamic Smogorinsky model (SDSM) 



for large-eddy simulation in a series of articles 
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191 ]. They construct a function of Smogorinsky 



model coefficients at the filter and test-filter scales, 



/(A) = C sm (aA)/C sm (A), 
where the expression of / (A) can be get through priori knowledge. Thus ([7]) can be rewrite as 



(8) 



M ij =/(A)(QA) 2 |5|5' ii -A i | ( S|^ 



IJ ! 



(9) 



Taking ([9]) into © , we can get the Smogorinsky model coefficient in the dissipation range. 

In this paper, we first propose a general constraint principle to construct the relation between 
two different scales. Then we give two concrete constraints through energ y spect rum and priori 
knowledge. It has been developed several types of energy spectrum so far 
we choose Kovasznay spectrum 



20h23I]. In this paper 



mixed models 



241 ] . and apply the two constrained condition to single model and 



2614281] . The results of LES with any one of the two constraints show some superiority 



than other models. 



II. THE GENERAL CONSTRAINT 



As is well known that traditional LES is done in the inertial range, and the SGS dissipation 



e equals to total dissipation approximately. In general, total dissipation is made up of two 



parts, the SGS dissipation and the viscous dissipation e v at the filter scale of LES, 

-SGS 



£ 



(10) 



And we have an identical equation on SGS dissipation at scale A and a A, 



(11) 



e^ GS £(aA) SGS ' 
Then we can construct the relation of e s and e v , 

£ SGS = p^uy ( 12 ) 

Thus, putting (|12p into (jlip . we can get the identical equation as following, 

(uj Sjj ) _ ( Tjj Sij ) 



Uj cxnu j.^. „u replace and T^, 



(13) 



Taking model stress and T i " lod to replace r« and Ty, we can get 



(14) 



Equation (|14p is the general principle of constraint. Next, we will give two concrete form of F {e u ) 
based on the existing spectrum theory and priori data. 

A. Constrained condition based on spectrum theory 

Kovasznay (1948) proposed a energy spectrum, and it supports a energy distribution with wave 
number from inertial range to dissipation range. The form of Kovasznay spectrum is 

4/3 N 



m 



^H(i) ' * = £ )1/4 - (15) 



The viscous dissipation at filter scale A can be expressed as 



fKc 

e\ = 2v \ k?E(k)dk, k c = ir/A. (16) 
J o 



From (1151) and (1161) we know that 



*=^(i-(i-£(£f)')- (17) 

At filter scale A, (jiop can be written as 

e = e £jGS + e £. (18) 

Thus, 

-v v i C SGS 

. *7A = (19) 



£ £^+4 GS 1 + 7?a' £a 
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From ([15]) . ([I7D and flE]), we can get 



4/3 N 



where, 



k r n f v 



k d A U 



3\ 1/4 



2 vv 



7T ( V 2 ■ 2v(S i jS lJ/ 



1/2 



A ( < ^(SijSij l 
TTU 1 ' 2 



A (2(SijSij) 



1/4 



1/4 



1/4 



i, \ 1/4 



-1/4 



Putting (f2Tj) into (p0|) . the result can be got 

1 



1 + ??A 



1 



/3a 



1-/3a(1 + ??a)" 1/3 



A 4 /3 



^2{SijSij) 



1/3- 



From (J22J), we know 



1. 



Therefore, the expression of F (e u ) can be written as 

F(e A ) = rj A e A . 



In general case, the equation can be got approximately that 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



where, denotes filtering at A and at scale aA. Thus, the a concrete form (|14|) can be expressed 
as 



(26) 



B. Constrained condition based on priori data fitting 



As described by Frisch et al. (1978), the (3 model quantifies the course of cascade 



2J. Similar 



with this theory, we can propose a supposition on energy dissipation of cascade. On some scale, the 
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viscous dissipation e v denotes the loss of kinetic energy and the SGS dissipation e SGS denotes the 
kinetic energy transferring into the next scale. Based on the scaling similarity of eddy structure, 
£ a/ £ a GS can ' 3e expressed by the function of A (A is non-dimensional number of filter scale) 

£a/4 GS = /(A). (27) 

As A — > 0, we know that e^ GS — > and e\ — > c (constant). While A is in the inertial range even 
the energy-containing range, e\ is ignorable in comparison with e^ GS . Therefore, we can give the 
expression of / (A), and two-term polynomial is used here, 



/(A) = C,(AT 1 + A^), 



(28) 



where, C v is the coefficient, and C u is variable with different kinetic viscous v. In order to check the 
hypothesis (j'28|) . we can give the result of priori as Fig. [JjFrom Fig. [TJ we can see that C u e^/e^ s 
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FIG. 1: (Color online) In the case of different v, the distribution of C v e^/e^" s with changing scale for priori 
Line with circle: v = 0.0006; Line with square: v = 0.001; Line with delta: v = 0.002; Line with RTringle: 
v = 0.005. 



can be expressed by an unified function. Thus, t proves that the hypothesis (|28p is reasonable. 
Then we choose one of the viscous v and fit the priori result of e\ / e^ GS . Thus we find the concrete 
form of 7 



7 1 



-2, 7 1 = -1. 



The constrained condition can also be expressed as below 

/Tmod c 



/ mod c. .\ 
\ ' ij °IJ I 



(29) 



(30) 



4/ (A-* + A-i) e^/((aA)-a + (aA)"i) ' 
III Al and III Bl introduce two concrete forms of (|14p . and then we will introduce a new mix dynamic 
model used in this paper. 
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III. THE CONSTRAINED LES MODEL 



In this article, we first construct two dynamic Smogorinsky models with constrained condition 
(|26p or (|30p (CDSM1, CDSM2) contrast with some single model as dynamic Smogorinsky model 
and SDSM. 

For model CDSM1, we can know that 

r™ od = C mod A 2 |S|>%, (31) 
T$° d = C teat (aA) 2 \S ij \S ij . (32) 



And the constrained condition is 

(CmodA 2 | , 5|'5«j'S'i_.) ^ (C tcst (aA) 2 \S\SijSij) 

Thus, we can get the expression of C moc j, 

(MjjLij) 



(33) 



where, 



Cmod " JM-M-y (34) 



M tJ = (|5 l^^ )??aA£ ° A A 2 |g|^ - A 2 \S\S tJ , (35) 



For model CDSM2, 



M .. = jlgl^XA^ + A- 1 )^ a2| ~~ 
(|5|^^-)(aA)-2 + (aA)-i) e ^ 

-A 2 |S|% (36) 

Different from other previous methods, the scale-dependent constrained condition can also ap- 
ply to mix models. Here we propose a constrained dynamic mixed helical SGS models. There are 
several SGS stress model commonly used in large eddy simulation of turbulence. The Smogorinsky 
model is a kind of eddy-viscosity model, and it can predict energy dissipation well, and has good 
numerical robustness and stability. But it can not predict the energy backscatter, and the corre- 
lation with the real SGS stress is not good . The nonlinear model, also called gradient or tensor 
eddy viscosity model, can successfully predict the structure of real stress and energy backscatter 
and does not need second filtering. Have thought of the advantages and disadvantages of the two 
models, we can construct a two-term mixed model T™ d • T™ od and T™ od can be expressed below, 

T ij° A = C ll/l« + °12h,iji ( 37 ) 



— C2lFi^j + C22F2,ij, (38) 

where f 1M = A 2 |S|^-, f 2>ij = A 2 g|; = (aA) 2 \S\Sij, ^ij = («A) 2 gg. Then, with 
the constrained condition (j26f) or (130j) , using (|3ip , we can construct the constrained dynamic mixed 
models CDMM1 and CDMM2, respectively. 



IV. THE ANALYSIS OF THE RESULTS 



Now we will give the analysis of the results of the forced steady and decaying homogeneous 
isotropic turbulence. In order to validate our model, a DNS of three-dimensional incompressible 
homogeneous isotropic turbulence is introduced here. It solves the forced N-S equations using a 
pseudo spectral code in a cubic box with periodic boundary conditions, and the numerical resolution 
is 512 3 . A Guassian random field is the initial flow condition, and it has an energy spectrum as, 



E (k) = Ak'U£ko b e k o , (39) 
where ko = 4.5786 and Uq = 0.715. The whole system is maintained by a constant energy input 



rate e = 0.1 in the first two wave number shells 



29]. 



To validate the constrained condition, we take four types of Smogorinsky models 
(CDSM1,CDSM2,SDSM and DSM) for comparison firstly. In this paper the LES and DNS of 
forced turbulence are all performed starting from a Gaussian random field with an inertial energy 
spectrum as (|33p . and the kinematic viscosity is v = 0.001. In Fig. [21 we present the energy spectra 
of the four type single models for forced steady isotropic turbulence at different filter scales, (a) 
filter scale A = 37r/64 (b) filter scale A = 37r/128. The filter-scale of LES in Fig. [2]^a) is close 
to the inertial range and (b) close to the dissipation range. We can see that the energy spectra 
of CDSM2 is closer to the DNS spectra than other spectra at most of the wave-number in (a). 
The energy spectra of CDSM1 and CDSM2 shows better than SDSM and DSM particular in high 
wave- number, and SDSM is a little better than DSM in (b). 

Fig. [3] shows us the evolution of model coefficients (CDSM1, CDSM2, SDSM and DSM) with 
time in steady turbulence, (a) A = 37r/64 (b) A = 37r/128. In Fig. [3l we can see that the 
distribution of the four model coefficients varied drastically with different filter scale. Near the 
inertial range (a), the model coefficients of CDSM1, SDSM and DSM are close, and near the 
dissipation range(b), the model coefficients of CDSM1 and CDSM2 are much closer. Totally, the 
value of DSM coefficient is smallest, and CDSM2 is biggest. 
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FIG. 2: (Color online) Energy spectra and for steady forced isotropic turbulence (a posteriori). Solid line: 
DNS; dashed line: CDSM1; DashDD line: CDSM2; line with diamond: SDSM; line with delta: DSM. (a) 
filter scale A = 3tt/64 (b) filter scale A = 3tt/128. 
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FIG. 3: (Color online) model coefficient evolution with time for steady forced isotropic turbulence for a 
posteriori, (a) at filter-scale A = 37r/64, (b) at filter-scale A = 37r/128. Line with square: CDSMI; Line 
with delta: CDSM2; Line with diamond: SDSM; line with circle: DSM. 



Shown in Fig. 0]are the dynamically calculated model and test model coefficients of CDSMI, 
(a) A = 37r/64 (b) A = 37r/128. In Fig. IU we can see that C mo d and C t e S t vary no more than 10% 
in (a), and while vary about 50%. From these results, we can see that taking C mo d and Ct es t the 
same value is not proper out of inertial range. 

Fig. [5] and Fig. [6] show us decay energy spectra of the four model evolution with time (t = 0, 6ro 
and 12to), starting from a fully developed statistical steady state, (a) A = 37r/64 (b) A = 37r/128. 
Here To is the initial large eddy turnover time scale, and is denned as tq = ir/ (kollo). The solid 
lines are DNS spectra, and dashed lines are energy spectra for the four models. From the two 
figures, we can see that CDSMI and CDSM2 display the similar trends and predict the energy 
spectra very well. And DSM model displays a worse status than other models. 

Through the analysis above all, we can say that single model with our constraints (CDSMI and 
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FIG. 4: (Color online) model and test coefficients of CDSMI distribution with time for steady forced isotropic 
turbulence for a posteriori, (a) at filter-scale A = 37r/64, (b) at filter-scale A = 37r/I28. Line with square: 
model coefficient; Line with delta: test coefficient. 




FIG. 5: (Color online) Energy spectra (filter-scale A = 37r/64) for decaying isotropic turbulence (a poste- 
riori), at t = 0, 6to and f2ro, where tq is the initial large eddy turnover time scale. Bold solid line: DNS; 
Dashed line: (a) CDSMf , (b) CDSM2, (c) SDSM, and (d) DSM. 

CDSM2) can improve the LES results of isotropic turbulence, generally. 

Then, let's see the case of mixed model. In Fig. [7J we present the energy spectra of three type 
mixed models (CDMM1, CDMM2, DMM) and DSM model at different filter scale, (a) A = 3vr/64 
(b) A = 3vr/128. In (a), the spectrum of CDMM1 shows closer to the DNS spectrum, and CDSMI 
and DSM show the similar trend, and are better than DMM in most wave-number. In (b) CDMM1 
and CDMM2 predict the energy spectra better than other models particularly in high wave-number. 
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FIG. 6: (Color online) Energy spectra (filter-scale A = 37r/128) for decaying isotropic turbulence (a poste- 
riori), at t = 0, 6to and 12to, where To is the initial large eddy turnover time scale. Bold solid line: DNS; 
Dashed line: (a) CDSMI, (b) CDSM2, (c) SDSM, and (d) DSM. 




FIG. 7: (Color online) Energy spectra and for steady forced isotropic turbulence (a posteriori). Solid line: 
DNS; dashed line: CDMM1; DashDD line: CDMM2; line with diamond: DMM; line with delta: DSM. (a) 
filter scale A = 3tt/64 (b) filter scale A = 3tt/I28. 

In Fig. [8j we show the distribution of SGS dissipation with 6/rj for the four LES models 
(CDMM1 ,CDMM2, DMM and DSM) and DNS, where 6 is the consecutive filter scale and rj is 
the Kolmogrov scale. From Fig. [HJ we can see that the SGS energy dissipation of CDMM1 and 
CDMM2 are much closer to the DNS' than that of DMM and DSM in most of the range. 

In Fig. we display the probability density functions (PDF) of the four LES models at the 
filter scale A = 3-7r/64, and the PDF of DNS is taken as comparison. As the limitation of DSM 
model, it can not predict the backscatters. Here, we can see that CDMM1, CDMM2 and DMM 
can capture the backscatters. Particulallywe also can be easy to see the PDFs of CDMM1 and 
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FIG. 8: (Color online) The SGS energy dissipation versus S/r] (a priori). Bold line: DNS; Line with square: 
CDMM1; Line with delta: CDMM2; Line with diamond: DMM; Line with circle: DSM. 
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FIG. 9: (Color online) Different probability density functions for the energy flux of DNS and other SGS 
models at the filter scale A = 3ir/64 (a posteriori) Bold line: DNS; Dashed line: CDMM1; Line with square: 
CDMM2; Line with delta: DMM; DashDD line: DSM. 

CDMM2 fit the PDF of DNS very well. Through a priori test, we can get the results that model 
stress weight rj^ od of CDMM1, CDMM2 and DMM have a high correlation with the real stress 
weight 7i2 over 85%, and the result of DSM has a low correlation with the real case not more than 
30%. 

We also display the distribution of high-order velocity increment with r/A in Fig. (TUJ where r 
is the separation distance and A is the filter scale A = 37r/64. From the figure (a) (b) (c), we can 
see that CDMM1, CDMM2 and DSM are all close to DNS very well with no obvious difference 
in the range r/A ^ 3. While in range of r/A sC 3, CDMM1, CDMM2 display better than DSM, 
but are still far from the result of DNS. In general, CDMM1, CDMM2 show better than DSM and 
DMM. 
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FIG. 10: (Color online) The distribution of high-order moments of longitudinal velocity increment with the 
separation distance r, where A is the filter scale, (a) Sf , (b) £g , (c) S$ . 

V. CONCLUSION 

In this paper, we first put forward a general LES method to simulate the turbulent flows. Using 
the method, we can ignore the assumption of scale invariance, and we can also do LES out of the 
inertial range. The main idea of the method is construct a relation of SGS dissipation and viscous 
dissipation between two scale, and we use such a relationship as an additional constraint condition 
in the dynamic models to do large eddy simulation. 

Through a priori and a posteriori tests, the single model CDSM1 and CDSM2 can predict 
energy spectra well in different range, they give the more proper model and test model coefficients. 
While the mixed model CDMM1 and CDMM2 also can predict energy spectra well in different 
range, and at the same time they predict the SGS energy dissipation more accurate, the PDF of 
energy flux much closer to the real case, and capturing the energy backscatter very well. Using the 
two models, the model stress has a higher correlation with the real SGS stress, and predicting the 
structure of the turbulence has a little progress. 

Contrast with other constrained LES method, our method gives up the character of scale in- 
variance and could also do LES well out of inertial range and even apply into simulating other type 
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turbulence, such as turbulence in the boundary lay, compressible turbulence, etc. 
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